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ABSTRACT 



Aims. The spatial distribution of the SMC stellar component is investigated from 2MASS data. The morphology of the different age 
populations is presented. The center of the distribution is calculated and compared with previous estimations. The rotation of the 
stellar content and possible consequence of dark matter presence are discussed. 

Methods. The different stellar populations are identified through a CMD diagram of the 2MASS data. Isopleth contour maps are 
produced in every case, to reveal the spatial distribution. The derived density profiles are discussed. 

Results. The older stellar population follows an exponential profile at projected diameters of about 5 kpc (~5") for the major axis 
and ~ 4 kpc for the minor axis, centred at RA: O'' 51""", Dec: -73° 7' (J2000.0). The centre coordinates are found the same for 
all the different age population maps and are in good accordance with the kinematical centre of the SMC. However they are found 
considerably different from the coordinates of the centre of the gas distribution. The fact that the older population found on an 
exponential disk, gives evidence that the stellar content is rotating, with a possible consequence of dark matter presence. The strong 
interactions between the MCs and the Milky Way might explain the difference in the distributions of the stellar and gas components. 
The lack in the observed velocity element, that implies absence of rotation, and contradicts with the consequences of exponential 
profile of the stellar component, may also be a result of the gravitational interactions. 
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1. Introduction 

SMC is classified as an irregular dwarf galaxy, characterized 
by a pronounced central feature of a Bar and an eastern ex- 
tension called the Wing. There has been a lot of investigation 
carried out studying the distribution of different stellar popula- 
tions of this galaxy. Initially, Freeman ( 1970) has calculated the 
scale length of the exponential component of the SMC from B- 
mag nitude photometry. He has fou nd a '=0.63 kpc. A long time 
later, iGardiner & HawkinsI (1199 ll) have found for a northern di- 
rection an exponential profile with a ' = 1.2 kpc, describing well 
the distribution up to r=6kpc. However this value was calculated 
over the observed projected radial distribution assuming a spher- 
ica Uy symmetric halo structur e. 

[Morgan & Hatzidimitrioul ( Il995h have carried out a spectro- 
scopic survey of carbon stars in the outer parts of the SMC. 
Regarding the spatial distribution of these stars, they conclude 
that at most surface densities the SMC appears elliptical with 
major axis parallel to the Bar, but there is a significant northward 
distortion of the outer contours. The diameter of the outermost 
contour is ~10''-12''. The authors also underline the existence of 
a spiral-arm-like extension of the carbon stars distribution south- 
wards at (projected) distances larger than 4" (~4 kpc) from the 
optical centre, although it's origin is not clear 

Counts of sources towards the Magellanic C louds from 
DENIS near-infrared survey have been studied by ICioni et al] 
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(l2000h . Their investigation was based on the differentiation in 
the (I-J, I) Colour Magnitude Diagram (CMD) in order to distin- 
guish between three groups of objects with different mean ages. 
The spatial distribution of the three age groups is found quite 
different the youngest stars exhibit an irregular structure while 
the older stars are smoothly and regularly distributed. The AGB 
and RGB stars are characterised by a regular, but double peaked, 
structure, an offset from the HI distribution, and a mean age dif- 
ference of the two maxima. 

iMaragoudaki et all (|2001|) have studied the spatial distribu- 
tion of the SMC stellar population according to their age, based 
on optical data. They have shown that the older stellar popula- 
tion shows a rather regular and smooth distribution, while the 
youngest stellar component (age less than 8x10^ yr) appears 
mainly along the northeast-southwest direction forming the Bar 

IStanimirovic et al.1 (|2004|) in a detailed study of the HI kine- 
matics conclude that the HI velocity field shows a large velocity 
gradient from the south west to the north east. The isovelocity 
contours of this velocity field show some symmetry, suggestive 
of a differential rotation. Some large-scale distortions in this ve- 
locity field are easily visible but could be related to positions of 
several supergiant shells. In the same article the authors conclude 
that, in contrast to HI distribution the old stellar populations ap- 
pear to have a spheroidal spatial distribution and a total absence 
of rotation. 

In order to further investigate this discrepancy between the 
gas and stellar component distribution, we have studied the dis- 
tribution of the various stellar components of the SMC from near 
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Fig. 1. Colour Magnitude Diagram K vs. J-K, of the SMC, from 
2MASS 

IR surveys. We know that IR observations characterize the mass 
of the galaxy better than those in shorter wavelengths where ab- 
sorption masks a large part of the stellar component. A detailed 
study of the stellar component as it is revealed from 2MASS 
is given in section 2. A comparison with available radio data is 
given in section 3, implementing the same approach used for the 
near IR data. Discussion and conclusions are given in section 4. 

2. 2MASS data 

The near-infrared (NIR) images were obtained by the Two- 
Micron All-Sky Survey (2MASS), an ongoing effort to map 
the entire sky at J-band (1.25 microns), H-band (1.65 microns) 
and K-band (2.17 microns) wavelengths. The 2MASS survey is 
led by the University of Massachusetts, with all data and im- 
ages processed at Caltech's Infrared Processing and Analysis 
Centre (IPAC). The survey utilizes two nearly identical 1.3- 
meter telescopes located at Mount Hopkins (Arizona) and at 
Cerro Tololo (Chile). While the pixel size is 2.0 arcseconds, 
the survey strategy of over-sampling yields an effective reso- 
lution of about 1 arcsecond. Exposure times are 7.8 seconds. 
2MASS has uniformly scanned the entire sky in three near- 
infrared bands to detect and characterize point sources brighter 
than about 1 mJy in each band, with signal-to-noise ratio (SNR) 
greater than 10. This has achieved an 80,000-fold improve- 
ment in sensitivity relative to earlier surveys. The data used 
here are from the 2MASS All-Sky Point Source Catalog (PSC) 
(fp_psc), at IPAC Infrared Science Archive (IRSA), Caltech/JPL 
(jhttp://irsa. ipac.caltech.edu/applications/Gator). 

IRSA gives choices for the origin of the default magnitudes 
and uncertainties in each band. We have requested only sources 
with satisfying quality of photometry for all J, H and K band. 
In order to accomplish that we choose the case where: the de- 
fault magnitude is derived from a profile-fitting measurement 
made on the 1.3 sec "Read_2" exposures. The profile-fit mag- 
nitudes are normalized to curve-of-growth-corrected aperture 
magnitudes. This is the most common type in the PSC, and is 
used for sources that have no saturated pixels in any of the 1.3 
sec exposures. 

The examined area of SMC is 6.9''xl2". Data from a 
small near by field of 2.1"x 1° are also obtained in order 
to be used later on for calculating the Galactic contribution. 
The coordinates of this near by field are 47.5<R.A.<60 and 
-68<Dec.<-67. It was selected to be close, but out of the 
SMC field, in order to represent well the distribution of 



background stars near SMC. The dimensions ware chosen 
to be big enough to make statistics for the SMC field. The 
initial catalogues contain 359234 stars for the SMC area 
and 10444 stars for the small nearby field. Total photomet- 
ric uncertainty for 2MASS data is given by IRSA (details at 
http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec4_5e.h tmll l. 
Sources with total photometric uncertainty greater or equal 
to 0.2 mag for any of the J, H and K bands were excluded, 
providing more reliable photometric measurements. Finally the 
catalogues used contain 293330 stars for the SMC field and 
8428 stars for the small nearby field. 



2.1. Identifying Stellar populations 

Major stellar populations are identified based on matching fea- 
tures of the observed colour-magnitude diagram with expected 
positions of known popu lations. Based on Hiparchos data of 
SMC stars. iGavrasI (l2003l) has determined the criteria for the dif- 
ferent star types on the K versus J-K CMD of the SMC. The A, 
B type stars are found at -1 < J-K < 0.2 and K < 15. The F and 
G type have 0.2 < J-K < 0.45 and K < 15. The K, M and faint 
Carbon Stars meet the criteria of 0.45< J-K < 1.3 and K < 15 
and finally Carbon stars have J-K > 1.3 and K < 15. Applying 
the criteria above, to the 293330 stars for the SMC field we find 
2618 A, B type stai's, 30042 of F and G type, 177471 K, M and 
faint Carbon Stars and 4705 Carbon stars. 

Two CMDs: the K vs J-K for the SMC area and the nearby 
field were produced. Both diagrams are divided in a grid with 
the same number of cells (22x16). The cells have dimensions 
that provide a fine grid, while they include enough number of 
stars to make statistics. The grids had the same number of cells 
and same starting and ending points for each axis, in order to 
compare numbers of stars with same K and J-K values. The cells 
of the nearby field were normalized to the equal area of the cells 
of the SMC area and their stars considered as background stars 
were subtracted randomly. The CMD of K versus J-K for the 
SMC field after subtraction of th e background contribution is 
shown in Fig[T] The criteria set bv lGavrasI (120031) for J-K values 
are indicated by the dashed lines. 

Since the different populations were defined, star counts 
were performed in a grid of 288x500 pixels for the SMC region. 
The isopleths contour maps of the SMC field that were produced 
for each of the four cases can be seen in Fig|2] It is very inter- 
esting to point out how well the old population (bottom right in 
Fig© reveals the center of mass of this stellar population. 

2.2. Distribution of late type stars 

In the case of K, M and faint Carbon stars the distribution re- 
veals the elliptical structure of a disk. The number density of the 
stars versus the distance from the centre has been estimated. The 
isopleths contour map of the SMC field was reproduced (FigO 
with bigger steps between the contours in order to allow more 
accurate selection of the image points that will be taken under 
consideration. The outermost contour is the same as in Fig|2]and 
will be the one specifying the diameter of each structure. We se- 
lected 4 directions a, b, c and d along the "minor" and "major" 
axis to derive their density profiles. The "major" axis is the one 
following the direction of the Bar, from north-east to south-west, 
while the "minor" axis is directing from south-east to north-west. 

Fitting exponential curves to the diagrams of number 
density of sources (N) versus (projected) distance from the 
centre of the galaxy (r), with exponential functions was 
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Fig. 2. Isopleths contour maps of the SMC different stellar populations, from 2MASS. Top left: A and B type stars, top right: F and 
G, bottom left: Carbon stars and bottom right: K, M and faint Carbon. 



very satisfying. The best fit functions are Y=15.9±2.2*exp(- 
1.0+0.5*X)-0.8+3.1, Y=19.7+6.8*exp(-0.4+0.3*X)-4.7+8.1, 
Y=16.6+1.8*exp(-1.5+0.6*X)+0.3+2.0, Y=17.8+2.0*exp(- 
1.3±0.5*X)-0. 1+2.1, for the directions a, b, c and d respectively. 
The values of a"' are 1.0, 2.5, 0.7 and 0.8 derived for 
each direction respectively, corresponding to a mean value 
of a"'= 1.25 kpc. For the power-law scenario, the best fit 
functions are Y=5.0+0.5*X(-''-^±"'>, Y=6.8+0.5*X(-"<'±°i', 
Y=4.2+0.6*X(-°-^±" ", Y=5.0+0.5*X*-" '^±° i', for the directions 
a, b, c and d respectively. The chi-square parameters are 0.04, 
0.03, 0.09 and 0.11 in the case of exponential fitting curves 
for the directions a, b, c and d respectively. The correspond- 
ing values for the power-law case are 1.6, 2.4, 1.0 and 2.4, 
indicating that exponential fitting represents better the data. 
Consequently the mass distribution derived here is described by 
an exponential disk as illustrated in Fig. |4] with the best fitting 
curves. Summing up, 2MASS reveals an exponential profile for 
the older stellar population with a (projected) diameter ~ 3.8 
kpc for the "minor" axis and ~ 4.8 kpc for the "major" axis. 
The center of the elliptical distribution is best determined by the 
right bottom map of Fig. H] at RA: O'' 51'", Dec: -73" 7'.2, that 
represents the K, M and faint carbon stars. 

The exponential disk model for K, M and faint Carbon stars, 
has to face the fact that the values of a ' are not found the same 
for the semi-axis's, and the fact that exponential along direction 
(d) is measured to be "steeper" than the exponential along di- 
rection (a). The observed asymmetry might be explained by the 



inclination of the SMC, although deprojection procedure is not 
clear enough in the case of this galaxy. However further investi- 
gation should put light on this issue. 

3. Radio data 

A radio image of SMC at 1400 GHz (21cm) was obtained 
from Uli Klein (private communication ) based on t he dat a 
that were derived from the work by Havnes et al.l (Il986l) . 
The data represent well the extended radio envelope of the 
SMC, well beyond any site of obvious active star formation 
and are appropriate for investigating the gas distribution 
at the SMC. The radial profile of this image is presented 
in Fig. |5] Using the previous methodology we created di- 
agrams of flux density (F) versus distance from the centre 
(r), and fitted data with exponential and power-law curves 
in order to reveal the most efficient mathematical descrip- 
tion of this distribution. The best fit exponential functions 
we found are Y=1208.1+2.0*exp(-2.2+0.01*X)-H30.7+1.9, 
Y=17243.6±270.6*exp(-3.7+0.02*X)-i-28.2+2.1, 
Y=1801.9±21.2*exp(-1.33±0.02*X)-585.6+14.6, and 
Y=248641.0±6690.8*exp(-4.5±0.02*X)4-60.13+1.6, for 
the directions a, b, c and d respectively. The values of 
a"' are 0.45, 0.27, 0.75 and 0.22 derived for each direc- 
tion respectively, corresponding to a mean value of a"' = 
0.42 kpc. For the power-law case, the best fit functions 
are Y=241.6±0.5*X(-°-'*=^""°i\ Y=456.9+0.5*X(-3-4±ooi\ 
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Fig. 3. Isopleths contour map of the SMC K, M and faint Carbon Stars, from 2MASS. Lines indicate the selected directions. The 
"major" axis of the galaxy is the one following the direction of the Bar, from north-east to south-west, while the "minor" axis is 
directing from south-east to north-west. 



Y=191.7+0.5*X(-'^ ''±" "''2), Y=3087.4+10.8*X'-5-5±°°", for the 
directions a, b, c and d respectively. The chi-square parameters 
are 590.0, 390.8, 33.8 and 471.9 in the case of exponential 
fitting curves for the directions a, b, c and d respectively, 
whereas the corresponding values for the power-law case are 
5542.13, 11285.9, 188.0 and 419.5. The comparison of the two 
sets of parameters allow us to suggest that exponential fitting 
represents better the distribution. We notice however, that in for 
direction (d) power-law is more efficient the Consequently the 
mass distribution derived here is described by an exponential 
disk as illustrated in Fig. |6] with the best fitting curves. 

These led us to a mean value of a"'= 0.42 kpc for the radio 
data. The centre is found at RA: O'' 55" 7^2, Dec: -72" 47' 59". 
The HI exponential profile has a scale length of 0.45 kpc along 
the south-east direction (a), 0.27 kpc along the north-east direc- 
tion (b), 0.75 kpc along the to north-west direction (c) and 0.22 
kpc along the south-west direction (d). Thes e coordinates of the 
center are in accordance with the results of IStanimirovic et"an 
(12004 derived from higher resolution data, and confirms that 
the differences found here between the central coordinates of the 
stellar and the gas component are real. 



4. Discussion & Conclusions 

The distribution of the stellar components in the near IR 
(2MASS) has been used to investigate the radial distribution of 
the stellar mass in the SMC. The isopleths in near IR wave- 



lengths are expected to represent more accurately the distribu- 
tion of stars since the observations are less affected by absorp- 
tion. 

In the case of the 2MASS data we set the criteria for col- 
lecting the different stellar populations and present their spatial 
distributions. It has been found that exponential profile fits bet- 
ter the distribution of the older population of the SMC. In Fig. 
|2]the bottom right map of K, M and faint carbon stars, reveals 
extremely accurately the center of the mass distribution of the 
oldest stellar population of the SMC, the coordinates of this cen- 
ter are found to be the same for all cases, independently the age 
of population they represent. 

This centre of the isopleths is found to be RA: O'' 51""", Dec: 
-73° 7' and the mean value for scale length is a"' = 1.25 kpc. 



4. 1 . The center offset 

The important characteristic of the distribution of the older stel- 
lar population of the SMC is found to be the offset of the cen- 
ter of their distribution from the centre of the distribution found 
by HI data dStanimirovic et all |2004j). In a more recent work, 
iPiatek et all (l2007h have calculated the kinematical centre of the 
SMC at RA: O'' 52^.8, Dec: -72° 30' from HST data being in 
good agreement with our calculations. They also explain that 
the presence of a Bar or a strong tidal disturbance can infer ex- 
tra streaming motions to a system and thus old and young pop- 
ulations can have distinct kinematics. The strong interactions 
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Fig. 4. Number of stars per pixel (N) vs. distance from the centre (r) (star symbols), fitted with exponential curves (solid lines) and 
power-law curves (dashed lines), for 2MASS data for the a, b, c and d directions respectively. 



of the SMC with the LMC and the Milky\ yay that have prob- 
ably c ontributed to the formation of the Bar ('Marag oudaki et aU 
1200 ih could also have as a result the observed difference in the 
centres of the distributions of the stellar and the gas compo- 
nent as well as the small ra dial velocity gradient among SMC 
stars (lBekki&ChibaLl2008h . that contradicts with the exponen- 
tial profile of the stellar component. 

IStanimirovic et al] ( 12004 ^ suggest that according to SMCs 
rotation curve a dark matter halo is not needed to explain its dy- 
namics. However reconsidering the values of the center of the 
SMC, adopting the center of the mass derived above, when cal- 
culating rotational velocity profiles may be important consider- 
ing the issue of dark matter presence in this galaxy. SMC is a 
dwarf irregu lar galaxy, which are expected to be dark matter 
dominated ( Cote et al.|[2000t iBegum et all l2006t IStrigari et al] 
[2OO8). However, rotational velocities do not support this case 
for SMC, probably because they have been affected by the inter- 
actions of the Magellanic system with the Milky way. 

4.2. The scale length values 

As mentioned before, lFreemanl(ll970h has investigated the struc- 
ture of the SMC through bright stars distribution and has found 
that their elliptical component is characterised by a scale length 
of a '=0.63 kpc. We find an exponential profile for the older 
stellar population with a ' = 1.25 kpc, while the gas exponential 
disk has a '= 0.42 kpc. These three values represent well the 



structure of three different galactic components and are in accor- 
dance with the fact that gas (radio data) is more concentrated on 
a thinner disk, younger stars (best represented by B-magnitude 
distribution) lay on a thicker disk and finally oldest stars reach 
bigger distances from the center 
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Fig. 6. Flux density (F) versus distance from the centre (r) (star symbols), fitted with exponential curves (sohd lines) and power-law 
curves (dashed hnes), for radio data for the a, b, c and d directions respectively. 



